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j!  INTRODUCTION 

•"  'Smoke  screening  as  a  battlefield  countermeasure  is  highly 
dependent  upon  a  number  of  atmospheric  parameters,  not  a  function  of 
munition  expenditures  alone.  Consideration  must  be  given  to  ambient 
atmospheric  conditions,  forecasts,  and  the  aptly  named  "fog  of  war" 
and  its  affects  upon  the  optical  characteristics  of  the  atmosphere. 
Determination  of  the  optimum  smoke  density  to  render  the  atmosphere 
opaque  to  energy  from  the  visible  through  the  far  infrared  band  of 
the  electromagnetic  spectrum  can  be  accomplished  by  an  algorithm 
based  upon  atmospheric  optics  and  turbulent  diffusion  hypotheses. 
The  algorithm  KWIK  was  initially  conceived  as  a  munitions 
expenditures  model  but  rendered  versatile  enough  to  be  used  for  large 
area  screening  and  the  "seeability”  on  a  battlefield. 


The  philosophy  behind  the  design  and  implementation  of  the 
algorithm  was  to  keep  it  simple  by  using  conventional,  well-proven 
formulae  and  calculating  variables  within  the  program  from  regression 
fits  to  easily  obtained  or  measured  atmospheric  parameters.  The 
model  was  formatted  In  a  modular  sense  such  that  as  better  ecperi- 
mental  results  became  available  changes  could  be  made  with  minimum 
effort. 
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BACKGROUND 

The  decision  to  develop  a  multispectral  obscuration  model 
immediately  led  to  a  straw  man  outline  with  secondary  and  tertiary 
options  for  the  structure  of  the  algorithm  to  assure  flexibility. 
The  advantages  and  disadvantages  of  atmospheric  optic  and  diffusion 
postulates  that  could  be  used  were  weighed  against  the  meteorological 
observations  that  would  be  available  on  a  battlefield.  A  determina¬ 
tion  was  made  that  eight  available  surface  and  sensible  weather 
observations  were  basic  to  the  development  of  a  workable  obscuration 
scheme.  As  a  consequence,  the  atmospheric  optics  model  selected  was 
one  proposed  by  Downs  [1],  being  compatible  with  the  required 
meteorological  data.  The  optics  model  was  then  complemented  by  a 
diffusion  approach  advocated  by  Gifford  [2].  The  complementary 
models  require  a  methodology  for  classifying  atmospheric  stability. 
The  six  category  Pasquill  [3]  scheme  was  selected  on  the  basis  of 
compatibility  with  the  input  parameters.  The  relationships  between 
the  Pasquill  categories  and  observational  data  and  the  eight  meteoro¬ 
logical  inputs  to  the  model  are  given  in  tables  la.  and  lb. 

TABLE  la.  RELATION  OF  PASQUILL  CATEGORIES  TO  WEATHER  CONDITIONS 


A  -  Extremely  unstable 

D 

-  Neutral 

B  -  Moderately  unstable 

E 

-  Slightly  stable 

C  -  Slightly  unstable 

F 

-  Moderately  stable 

Nighttime  Conditions 

Sur face 

Thin  Overcast 

Windspeed  Daytime  Insolation 

or  _>  4/8  <_3/8 

(m/s)  Strong  Moderate  Slight 

Cloudiness  Cloudiness 
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TABLE  lb.  METEOROLOGICAL  INPUT  PARAMETERS  FOR  STABILITY, 
OPTICS,  AND  DIFFUSION  CALCULATIONS 


Ceiling  height  (feet) 

Temperature  (degrees  F)  1 

Cloud  cover  (percent) 

Dew  point  (degrees  F)  f 

Visibility  (miles) 

Wind  direction  (degrees)  j. 

Precipitation  (yes  or  no) 

Windspeed  (knots)  | 

The  Gaussian  plume  and  puff  hypothesis  chosen  to  model  h 

semlcontinuous  and  quasi-instantaneous  sources,  uses  longitudinal, 
lateral,  and  vertical  dispersion  parameters  that  are  stability 
dependent  and,  in  the  case  of  vertical  diffusion,  a  function  of  the 

aerodynamic  roughness  of  the  earth's  surface.  The  dispersion  i 

parameters  are  calculated  using  power  laws  attributed  to  Pasqulll  , 

[3] .  Aerodynamic  roughness  Is  evaluated  using  the  approach  of  Rung  1 

[4] ,  based  upon  the  height  of  the  roughness  elements,  l.e.,  trees,  r 

bushes,  grasses,  etc.  h 

The  straw  man  approach  highlighted  additional  factors  that  • 

would  Influence  obscuration  on  the  battlefield  including  (1)  relative 
hunidlty  effects  upon  the  hygroscopic  characteristics  of  the  smoke 
aerosols,  (2)  the  need  for  optimum  impact  separations  of  artillery-  ! 

delivered  smoke  projectiles,  (3)  rate  of  fire  calculations  for  effi¬ 
cient  dissemination  of  chemical  smokes,  and  (4)  munition  expenditure  i 

estimates  used  for  planning  future  operations.  | 

THE  FOG  OF  WAR  | 

i 

The  visual  range  in  the  atmosphere  is  directly  affected  by  l 

natural  aerosol  concentrations  and  particle  size  distributions  occur-  | 

ring  over  any  optical  path.  Visibility  may  be  reduced  by  dry  haze,  I 

wet  or  relative  humidity  haze,  fog,  or  air  pollution.  Battlefield  | 

visibility  can  be  compromised  by  the  additive  effects  of  the  fog  of 
war,  l.e.,  pollution  induced  by  deliberate  smoke  screening,  dust 
thrcwn  up  by  the  mass  movement  of  heavy  vehicles  and  intense 
artillery  barrages,  or  smokes  from  burning  vehicles.  The  reduction 
of  visibility  by  natural  aerosols  plus  the  fog  of  war  had  a  direct 
effect  upon  countermeasure  obscuration  used  to  deny  target  acquisi- 
t  ion. 

Battlefield  pollution  coupled  with  the  liquid  water  content  I 

of  the  atmosphere  control  the  aerosol  scattering  and  water  vapor  J 

absorption  of  visible  light  and  infrared  radiation  in  the  J 
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atmosphere.  The  attenuation  of  an  optical  path  by  scattering  and 
absorption  results  in  a  higher  threshold  level  of  detection.  The  net 
effect  is  a  need  for  less  smoke  to  render  an  optical  path  opaque. 
The  KWIK  algorithm  was  designed  to  take  these  phenomena  into  account, 
i.e. ,  aerosol  mass  concentrations  with  respect  to  battlefield 
visibllty  and  attenuation  by  water  vapor  as  a  function  of  relative 
humidity. 


The  inclusion  of  visibility  in  the  determination  of  extinc¬ 
tion  coefficients  for  finite  optical  paths  is  a  major  factor  in  the 
reduction  of  munitions  expenditures  necessary  to  establish  and  main¬ 
tain  obscuring  screens  on  a  battlefield.  Conversely,  as  discussed 
below,  the  KWIK  algorithm,  with  some  modifications,  can  be  used  to 
estimate  the  degradation  of  visibility  by  the  products  of  a  dirty 
battlefield. 

SMOKE  SCREENING  AND  OBSCURATION 

Countermeasure  obscuration  may  be  considered  from  two  view¬ 
points,  the  first  being  the  attenuation  of  near  horizontal  optical 
paths  where  the  crosswind  integrated  concentration  of  smokes  such  as 
white  phosphorous  or  zinc  chloride  is  of  prime  importance.  The 
second  is  large  area  screening  whereby  detection  of  ground  targets  by 
airbone  observers  or  remotely  piloted  vehicles  is  denied.  Here  the 
vertical  or  slant  path  Integrated  concentration  is  the  prime 
parameter.  The  determination  of  the  horizontal  or  vertical 
integrated  concentrations  necessary  to  obscure  an  optical  path  to  the 
threshold  level  is  almost  wholly  dependent  upon  the  basic  optical 
calculations.  The  optics  portion  of  KWIK  is  adapted  from  an  approach 
to  atmospheric  transmission  suggested  by  Downs  (lj.  Transmittance  of 
light  at  various  wavelengths  through  a  path  is  determined  by  calcu¬ 
lating  the  attenuation  due  to  absorption  by  water  vapor,  scattering 
by  haze  or  fog,  and  precipitation.  When  the  attenuation  due  to  atmo¬ 
spheric  conditions  is  known,  the  attenuation  due  to  smoke  that  is 
required  to  lower  tansmlttance  to  a  threshold  contrast  can  be 
computed. 


Absorption  is  attributable  to  the  amount  of  precipitable 
water  in  a  path,  assuming  the  water  vapor  concentration  in  the  atmo¬ 
sphere  Is  well  behaved  and  exhibits  a  scale  height  of  about  2  km. 
The  water  vapor  concentration  expressed  in  centimeters  pet  kilometer 
of  path  length  may  then  be  given  as: 

W  -  WQe"<L  sin  9 ,  (1) 
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where  WQ  is  the  precipitable  water  along  a  path  L  and  Q  is  the  angle 
between  the  horizontal  and  the  height  of  a  target  above  or  below  an 
observer.  WQ  is  computed  from  a  regression  equation  relating 

precipitable  water  and  dew  point  temperature  (Td): 

2  3 

wc  -  0.4477  +  0.0328Td  +  1.2(10)“3  Td  +  1.84(10)“5  Td  (2) 

Equation  (2)  was  fit  to  data  extracted  from  Downs  [1]  and  is 
considered  valid  for  any  meteorological  condition. 

The  amount  of  water  vapor  in  the  path,  W,  is  given  by: 

W  -  «  /L  8l"  e)/2  dL  .  (3) 

Transmission  through  the  absorbing  component  of  the  atmo¬ 
sphere  is  calculated  by  using  an  error  function  absorption  law  devel¬ 
oped  by  Elsasser  [5] 


T  -  1  -  erf  (z) 


(4) 


erf  (z) 


(5) 


where  z  =■  0.5  3  /iiw  ,  and  6  is  the  error  function  absorption 

coefficient  as  a  function  of  wavelength. 

Downs  (11  states  that  the  Elasasser  approach  is  unable  to 
correctly  address  long  wavelengths  and  suggests  using  an  approach 
described  by  Fisher  [6]  for  the  far  infrared  wavelengths,  with  the 
computation  of  transmission  due  to  absorption  by  water  vapor  given 
by: 

T  -  e-0.0681W  #  (6) 


Reduction  in  transmittance  due  to  attenuation  by  haze  and 
fog  can  be  calculated  by  using  Mie  theory.  Downs  fl]  indicates  that 
the  Mie  scattering  coefficient  decreases  with  altitude  such  that  its 
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behavior  can  only  be  estimated.  The  following  expressions  for 
(Mle  scattering  coefficient)  are,  at  best,  an  approximation  to  the 
behavior  of  the  versus  altitude  relationship 


o^fe"L  sin  0/4,1  *,  vR  ±G(X)  , 


(7) 


%  " 

^2 


o^e1,  sin  e  In  (0.1/c^f)  «  VR  <  G(  X)  ,  0  <  L  sin  0  £1  km 
-  0.128e-L  sln  9/4,1  ;  VR  <  G(X)  ,  1  km  <_L  sin  9  <  » 


°M  *  %  (°M2)  *  (fl) 

where  Is  the  extinction  coefficient  determined  from  a  linear 

regression  as  a  function  of  visibility  and  wavelength,  based  upon 
Downs'  evaluation.  VR  is  visibility,  and  G(  X)  Is  the  scale  height  of 
cfy.  G( X)  is  not  constant;  rather  it  Is  a  function  of  altitude,  visi¬ 
bility,  and  wavelength.  Transmission  along  a  path  with  attenuation 
can  be  determined  by  the  equation 

T  -  e_  /L  °M<L)  dL  ,  (9) 

o 

after  substituting  a  value  for  cfy  according  to  equations  (7)  and 
(8).  If  precipitation  is  Indicated  (by  Input  parameter),  then  the 
value  for  transmittance  In  equation  (8)  is  set  to  one  and  a  calcula¬ 
tion  is  made  for  attenuation  by  precipitation  Instead. 

Reduced  transmittance  owing  to  attenuation  by  precipitation 
can  be  obtained  from 


T  -  e“  /  <J  (L)  dL  ,  (10) 

o  r 

where  L  is  the  path  length  and  is  an  attenuation  coefficient 
determined  from  a  regression  equation  as  a  function  of  visibility  and 
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wavelength.  The  total  transmittance  along  an  optical  path  Is  the 
product  of  the  partial  transmlttances 

Ttotal  “  VhfVa  *  <U> 

where 


Tg  -  Transmittance  due  to  attenuation  by  absorption 

Thf  ■  Transmittance  due  to  attenuation  by  haze  and  fog 

Tp  -  Transmittance  due  to  attenuation  by  precipitation 

Tg  *  Transmittance  due  to  attenuation  by  smoke 

Tg  is  then  calculated  from  equation  (11),  and  the  desired  threshold 
contrast  of  transmittance  for  a  particular  wavelength  can  be  deter¬ 
mined  from: 


T 

s 


tc 


T  T..T 
a  hf  p 


(12) 


where  Ttc,  the  threshold  contrast,  ia  baaed  upon  the  Koschmieder  f7) 
theory  but  set  equal  to  0.10  for  visible  wavelengths  and  0.05  for 
infrared. 


The  approach  used  to  determine  the  line  of  sight  integrated 
concentration,  CL,  of  the  smoke  screen  necessary  to  attenuate  an 
optical  path  to  a  threshold  level  is  based  on  the  transmittance  of 
equation  (12)  as  a  function  of  CL.  The  CL-value  necessary  to 
attenuate  an  optical  path  to  deny  target  acquisition  may  be 
determined  from  the  Bouguer-Beer  law  written  as 


T 

s 


T  -aCL 

T"“  e 

o 


(11) 


where  I  is  the  illuminance  at  a  target,  I  is  the  illuminance  at  the 
light  source,  a  the  extinction  coefficient,  and  Tg  the  reciprocal  of 
attenuation.  Rearranging  and  solving  equation  (11)  yields 


CL  ■  If. 

-a 


(14) 


where  CL  is  the  minimum  integrated  smoke  cone  itration 
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Equation  (14)  Is  applicable  to  both  visible  and  Infrared 
wavelengths,  If  values  of  a  are  known.  Listed  In  table  2  are  average 
extinction  coefficients  for  the  two  screening  chemical  smokes  consid¬ 
ered  by  the  KWIK  algorithm,  l.e. ,  bulk  white  phosphorous  (WP)  and 
hexachloroethane  (HC)  (zinc  chloride)  for  the  Indicated  spectral 
bands. 


TABLE  2.  EXTINCTION  COEFFICIENTS  FOR  WHITE  PHOSPHOROUS 
AND  HEXACHLOROETHANE  GENERATED  SMOKES 


Spectral 

Band 

Wavelength 
(  um) 

Zinc  Chloride 

— j — n - 

m  g 

White  Phosphorous 

Visible 

0.4  to  0.7 

3.30 

2.46 

Near  IR 

0.75  to  1.2 

1.50 

1.50 

Mid  IR 

3  to  5 

0.12 

0.21 

Far  IR 

8  to  14 

0.05 

0.28 

The  development  of  diffusion  formulae  for  predicting  the 
obscuring  power  of  chemical  smokes  starts  with  the  assumption  of  an 
Instantaneous  point  source  of  material  diffusing  in  three  dimen¬ 
sions.  For  a  Gaussian  distribution  of  diffusion  taking  place  Inde¬ 
pendently  in  the  three  coordinate  directions,  the  equation  can  be 
stated  as 


X(x,y,z) 


(2  IT) 


-3/2  Q_ 


a  a  a 
x  y  z 


exp 


(15) 


where  x  Is  concentration  In  g  ra  J,  Oy  the  total  release  of  material 
In  g,  x,  y,  z  the  coordinate  directions,  V  the  mean  windspeed  In 
m  s~  ,  t  Is  time,  and  ox,  <jy,  az  the  dispersion  parameters. 

Integration  of  equation  (15)  yields  the  continuous  source 
equation  for  sources  and  receptors  near  the  ground  as 


it  V  a  a 
y  z 


exp 


[-*(*)’•(*)’] 


(16) 
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where  Q  is  a  time  rate  of  release.  If  equation  (16)  is  integrated, 
the  result  is  the  crosswind  integrated  concentration  (CWIC)  of  the 
plume 

'  77  '^["Kt)2]  (17) 

which  is  the  basic  form  for  estimating  obscuration  from  sources 
generated  by  the  HC  smoke  mix.  The  vertical  dispersion  parameter  oz 
is  determined  using  the  Pasquill  [3]  power  law  in  the  form 

o2  -  cxd  .  (18) 

The  coefficient  and  index  values  will  be  discussed  at  the  end  of  the 
section. 


Screening  and  obscuration  are  not  restricted  to  the  along 
wind  case  of  equation  (16);  consequently,  a  wind  direction  correction 
factor  must  be  considered  for  head,  tail,  and  quartering  wind  condi¬ 
tions.  The  correction  factor  may  be  derived  from  considerations 
concerned  with  finite  line  sources,  cumulative  effects  of  multiple 
sources,  and  discrete  point  line  sources.  The  correction  factor  is 
eliptical  and  given  by 


2  2 
m  n 


2  2  2  2 
m  sin  9  +  n  cos  9 


(19) 


where  m  «  3.71  and  n  »  1,  the  eliptical  semi-axes,  and  9  the  angle 
between  the  mean  wind  direction  and  the  optical  path. 

The  exponential  term  on  the  right  hand  side  of  equation 
(17)  is  only  partially  sensitive  to  stability  and  downwind  travel 
distances  associated  with  smoke  screening.  Consequently,  numerical 
evaluation  for  the  six  stability  categories  and  along  wind  distances 
of  50  to  150  m  shows  that  exp[-l/2(z/oz)  ]  may  be  taken  as  constant 
and  set  equal  to  0.916. 

The  assumption  can  be  made  that  the  line  of  sight 
Integrated  concentration  CL  of  the  screening  aerosol  calculated  from 
equation  (14)  is  equal  to  the  value  of  Xcyjc  evaluated  from  equation 
(17).  This  allows  equation  (17)  to  be  rearranged,  after  substitution 
of  equation  (18)  and  consideration  of  the  ramifications  of  equation 
(19)  and  solved  for  the  along  wind  travel  distance  x.  Included  in 
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the  solution  are  the  source  efficiency  term  X  and  the  relative  hunid- 
ity-related  yield  factor  fi,  which  serve  to  modify  the  source  strength 
Q.  Thus 


X 


-1 1  0.731  X  ft  Q  | 
L  C  V  CL  J 


.-1 


(20) 


after  combining  all  the  constants.  The  integrated  concentration 
required  for  obscuration  reaches  a  minimum  at  the  distance  x  downwind 
from  the  source.  Accordingly,  x  is  the  calculated  impact  separation 
of  the  smoke  projectiles. 

Obscuration  calculations  for  quasi-instantaneous  sources, 
i.e.,  as  generated  by  bulk  WP  munitions  requires  that  the  integrated 
concentration  equation  be  written  in  the  form 


XCWIC 


X  0  Q„ 
_ _T 

wo  o 


where  the  dispersion  parameters 


and 


o-  are  not  to  be  confused 
ZI 


with  those  associated  with  a  continuous  source.  The  term  Z  repre¬ 
sents  the  height  of  the  puff  centroid  above  the  surface.  Owing  to 
the  nature  of  the  exponential  on  the  right  hand  side  of  equation  (21) 
which  requires  knowledge  of  two  downwind  travel  distances  and  three 
heights  for  solution,  the  integrated  concentrations  were  determined 
for  a  unit  downwind  distance  of  100  m. 

Bulk  WP  is  an  exothermal  smoke  source,  with  only  a  fraction 
of  the  total  material  available  for  screening.  The  major  portion  of 
the  phosphorous  smoke  is  transported  vertically  in  the  thermal 
plume.  An  analysis  of  available  data  yielded  the  efficiencies  or 
percent  of  payload  available  shown  in  table  3.  Also  tabulated  in 
table  3  are  numerical  evaluations  of  the  exponential  term  of  equation 
(21)  for  each  stability  category  (K  factors). 
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TABLE  3.  EFFICIENCIES  AND  K  FACTORS  FOR  BULK  WHITE  PHOSPHOROUS  SMOKE 
MUNITIONS  AS  A  FUNCTION  OF  PASQUILL  STABILITY  CATEGORY 


Pasquill 

Category 

K  Factor 

Ef  f iclency 

A 

0.4633 

0.07 

B 

0.3631 

0.10 

C 

0.2056 

0.14 

D 

0.0647 

0.28 

E  and  F 

0.0725 

555  0.  30 

This  allows  equation  (21)  to  be  restated  as 

k  x  a  0T 

*CWIC  =  T~a  o  * 

XI  ZI 

Projectile  impact  separations  can  now  be  determined  from 


x  =  100 


^WIC 

CL 


(22) 


(23) 


where  the  constant  is  the  unit  screen  length  and  CL  Is  calculated 
from  equation  (14). 


Large  area  screening  can  be  treated  by  integration  of  equa¬ 
tion  (16)  in  the  vertical  yielding 


*VIC 


(24) 


which  is  applicable  to  the  use  of  fog  oil  generators.  The  rendering 
opaque  of  large  areas  such  as  air  fields  is  dependent  upon  the  inher¬ 
ent  and  apparent  contrast  between  two  objects  on  the  ground  with 
respect  to  the  sky-ground  ratio  and  the  threshold  contrast  of  the 
objects.  The  problem  of  vision  looking  downward  from  aircraft  has 


or  A 


¥r'  "riVifirii'lirimnrtii  m  i 
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been  discussed  by  Duntley  [8]  and  Middleton  [9],  Middleton  suggests 
that  the  threshold  contrast  over  a  slant 

range  R  is  given  by 


e  =  C 

o 


/  B”  \  / 

'  bo  '  V 


-.-1 


exp  3.912  R/V 


(25) 


where  CQ  is  the  inherent  contrast,  Bm/B0  the  sky-ground  ratio,  and  Vr 
the  visual  range.  Values  of  e  so  determined  may  be  used  in  equation 
(14)  to  evaluate  slant  range  integrated  concentrations. 
Sherwood  [10]  found  that  ^*33  g  m~^  were  necessary  to 
screen  a  large  area  from  aerial  observation. 

The  obvious  use  of  equations  (24)  and  (25)  is  the  determi¬ 
nation  of  the  number  and  separation  of  oil  fog  generators  required  to 
screen  a  large  area.  Equation  (24)  can  be  manipulated  to  yield  the 
separation  distance  yg  as 


2  c 

In  (- 

l  \  0 

y 

\ 

V  V  o  X  _ 

y  VIC  J 

(26) 


and  the  number  of  generators  Ng  required  by  dividing  the  screen  width 
L  by  ys  plus  one,  or 
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(27) 


The  dispersion  coefficients  ay>  oz,  ox^,  and  oz^  are  based 
upon  the  Pasquill  [3]  power  laws  and  given  by 

°y  =  axb  (28) 

az  *  cxd  (29) 
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(J  =  a  +  0.667  cx^  (31) 

ZI  zo 

where  ax  and  az  are  the  Initial  dispersion  or  burst  functions. 
Downwind  °dispersl8n  parameters  for  the  quasl-lnstantaneous  sources 
are  approximately  two-thirds  of  those  for  continuous  sources  as  shown 
by  Pasquill  [3],  which  Is  reflected  by  equations  (30)  and  (31).  The 
coefficients  and  indices  for  the  dispersion  parameters  as  a  function 
of  stability  and  aerodynamic  roughness  are  tabulated  in  table  4  for 
three  roughness  lengths. 

TABLE  4.  COEFFICIENT  AND  INDEX  VALUES  FOR  THE 
POWER  LAW  DISPERSION  PARAMETERS 


Stability 

Category 

a 

b 

C 

zo  =  1 

d 

cm 

c  d 

zQ  =  10  cm 

c  d 

z0  =  100  era 

A 

0.40 

0.90 

0.154 

0.94 

0.279 

0.90 

0.615 

0.83 

B 

0.32 

0.90 

0.133 

0.89 

0.225 

0.85 

0.539 

0.77 

C 

0.22 

0.90 

0.121 

0.85 

0.213 

0.81 

0.533 

0.72 

D 

0.143 

0.90 

0.108 

0.81 

0.195 

0.76 

0.456 

0.68 

E 

0.102 

0.90 

0.078 

0.78 

0.139 

0.73 

0.348 

0.65 

F 

0.076 

0.90 

0.062 

0.72 

0.117 

0.67 

0.309 

0.58 

SEEABILITY  ON  THE  BATTLEFIELD 

A  battlefield  may  be  considered  as  being  mesometeorological 
in  scale,  i.e.,  areas  ranging  from  hundreds  to  thousands  of  square 
kilometers.  If  the  density  of  meteorological  observations  is  large 
and  timely,  the  optics  portion  of  KWIK  may  be  used  to  calculate  the 
attenuation  of  optical  paths  for  each  weather  observational  site. 
The  attenuations  may  be  plotted  and  analyzed  much  like  synoptic  data 
to  prepare  "seeabllity"  charts  for  a  battlefield.  Battle  plans  for 
future  engagements  with  estimates  of  munition  expenditures 
anticipated  number  of  burning  vehicles  and  vehicular  dust  conditions 
can  be  used  to  predict  visibility  conditions  which  may  be  used  for 
attenuation  forecasts.  Seeabllity,  prognostications  can  be  used  for 
planning  purposes,  i.e.,  what  weapons  system  will  be  effective  on  the 
next  day's  predicted  dirty  battlefield. 
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DISCUSSION 

The  complete  KWIK  smoke  obscuration  model  has  been 
discussed  by  Umstead,  Pena,  and  Hansen  [111,  including  the  develop¬ 
ment  of  the  scheme  to  establish  rates  of  fire,  impact  separations  in 
adverse  wind  conditions,  and  the  special  considerations  for  munitions 
expenditures  in  the  mid  and  far  infrared  regions  of  the  spectrum. 
Owing  to  space  limitations,  these  subjects  are  beyond  the  scope  of 
this  paper. 

The  basic  scheme  developed  for  artillery  delivered  smoke 
projectiles  assumes  that  smoke  obscuration  operations  will  be 
conducted  with  batteries  or  battalions  firing  in  an  open  sheaf  rather 
than  parallel  or  normal  sheaf  patterns.  This  concept,  coupled  with 
relative  hunidity  dependent  yield  factors  in  the  smoke  model,  will 
result  in  a  more  efficient  use  of  smoke  on  a  battlefield. 
Preliminary  studies  suggest  savings  in  munitions  up  to  20  to  30  per¬ 
cent  over  current  methods. 

The  use  of  the  KWIK  algorithm  for  large  area  screening 
operations  and  predicting  attenuation  degradation  on  the  dirty 
battlefield  does  not  detract  from  the  original  intent  of  the  model, 
but  enhances  its  capabilities.  The  outputs  of  these  offshoots  of  the 
primary  model  can  be  utilized  to  improve  the  munition  expenditure 
estimates  generated  for  countermeasure  obscuration  purposes.  Large 
area  screening  systems  utilizing  oil  fogs  are  only  usable  in  the 
visible  portion  of  the  spectrum.  Consequently,  the  portion  of  KWIK 
expressed  by  equations  (24),  (25),  and  (26)  only  apply  to  the  0.4-  to 
0.7 -vim  band. 

Owing  to  the  extinction  character istics  of  HC  smoke,  the 
approach  personified  by  equation  (20)  is  valid  only  in  the  visual  and 
near  infrared  portions  of  the  electromagnetic  spectrum.  The  superior 
characteristics  of  WP  allow  it  to  be  utilized  to  countermeasure 
devices  operating  in  the  mid  and  far  infrared  regions. 

CONCLUSIONS 

The  KWIK  obscuration  model  is  highly  versatile  and  has  been 
programmed  to  operate  on  a  variety  of  machines  ranging  from  program¬ 
mable  desk  calculators  to  digital  computers.  Three  versions  ate  in 
existence:  the  Fortran,  the  real-time  with  three  options,  and  a 
deferred-time  version  for  generating  munition  expenditure  estimate 
tables  based  upon  climatological  input  data.  The  algorithm  may  be 
used  in  threat  analysis  studies  or  as  a  subroutine  in  force  on  force 
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scenarios.  The  modular  concept  used  to  develop  RWIK  allows  greater 
flexibility  than  found  In  predecessor  smoke  models. 
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